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A bilayer system of two-dimensional electron gases in a perpendicular magnetic field exhibits rich 
phenomena. At total filling factor i/tot = 1, as one increases the layer separation, the bilayer system 
goes from an interlayer coherent exciton condensed state to an incoherent phase of, most likely, 
two decoupled composite-fermion Fermi liquids. Many questions still remain as to the nature of 
the transition between these two phases. Recent experiments have demonstrated that spin plays 
an important role in this transition. Assuming that there is a direct first order transition between 
the spin-polarized interlayer-coherent quantum Hall state and spin-partially-polarized composite 
Fermi liquid state, we calculate the phase boundary {d/l)c as a function of parallel magnetic field, 
NMR/heat pulse, temperature, and density imbalance, and compare with experimental results. 
Remarkably good agreement is found between theory and various experiments. 



I. INTRODUCTION 

A bilayer two dimensional electron gas in a strong per- 
pendicular magnetic field at total filling factor vtot = 1 
exhibits rich phenomena. An important tuning param- 
eter in this system is the ratio djl, where d is the ef- 
fective interlayer distance, and I is the magnetic length. 
At small djl, even with negligible tunneling, a remark- 
able bilayer quantum Hall state with interlayer phase co- 
herence emerges due to interlayer Coulomb interaction. 
This bilayer quantum Hall state can be described as a 
pseudospin ferromagnet where the layer index acts as 
the pseudospin, or an exciton condensate formed by in- 
terlayer particle-hole pairing^"'^. Many remarkable ex- 
perimental signatures of this phase predicted by theories 
have been observed in experiments, including enormous 
enhancement of zero bias interlayer tunneling"*, linearly 
dispersing Goldstone mode^, quantized Hall drag^, and 
vanishing resistance in counter flow^. However, there are 
still important discrepancies between theory and exper- 
iment. For example, the height of the interlayer tunnel- 
ing conductance is observed to be finite'*, while theories 
predict it to be infinite. Also, transport in counterflow 
experiments should be completely dissipationless under 
a critical temperature for phase coherence, but in experi- 
ments dissipationless counterflow is only seen in the zero- 
temperature limit''. The effect of quenched disorder is be- 
lieved to be crucial to reconcile these discrepancies^"^°, 
although a quantitative understanding is still lacking. 

The nature of the phase transition when djl is in- 
creased and the quantum Hall phase is destroyed is 
even more puzzling. In the limit djl — oo, each 
layer is at half-fllling, and they should behave as 
weakly-coupled composite Fermi liquids. Much progress 
have been made in understanding this phase using the 
Chern-Simons approach*^^^^ and the dipolar quasiparti- 
cle approach*^^* . Although we understand well both 
the coherent phase at djl — >■ and the composite 
Fermi liquid state at djl — oo, the transition between 
them has been shrouded in mystery. There have been 



many experimentaP^"'^^ and thcoreticaP^"''^ studies re- 
garding the nature of this transition. While some of 
these theoretical works point to a direct transition be- 
tween the two limiting phases, either continuous^^ or of 
flrst order'*^"'*'^, some other works predict the existence 
of various types of exotic intermediate phases, including 
translational symmetry broken phase^^"'^^'^^, composite 
fermion paired state^^'^°'''^, phase of coexisting compos- 
ite fcrmions and composite bosons^^'^®'"*^, and quantum 
disordered phases'**, etc. 

These theoretical works typically assume that the 
physical spin is fully polarized and hence irrelevant across 
the transition. However, recent experiments have shown 
that spin plays an important role in the transition. Ref. 
26 has found that by applying a NMR pulse or heat pulse 
to depolarize the nuclei and hence increasing the effective 
magnetic held coupled to the spin, the coherent phase is 
strengthened, and the phase boundary shifts to higher 
value of djl. Similar behavior has also been observed 
by applying a parallel magnetic field^''. These experi- 
mental results indicate that at least one of the phases 
involved in the transition is not fully polarized, and that 
the polarization changes significantly accross the transi- 
tion. The most likely possibility is that the incoherent 
composite Fermi liquid phase at large djl is only par- 
tially polarized, as shown by other experiments on single 
layer at = 1/2^°'^*. If the transition between the co- 
herent phase and the less 'polarized incoherent phase is a 
thermodynamic phase transition, it must be of first order: 
The magnetization is discontinuous across the transition, 
and, as the experiments of Ref. 26 found, the transition 
can be tuned using a Zeeman field which is conjugate to 
the magnetization. These two facts together imply the 
first order nature of the transition. An alternative to 
the thermodynamic transition scenario is a singularity- 
free quantum crossover as was suggested recently in Refs. 
47,48. 

In this work, we assume that the transition tuned by 
djl is a thermodynamic first-order transition between 
spin-polarized coherent vtot = 1 quantum Hall state and 
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partially-polarized composite Fermi liquid state, and de- 
rive the Clausius-Clapeyron relations for this system. 
The Clausius-Clapeyron relations will allow us to obtain 
the phase boundary shapes for the transition; a com- 
parison of these boundaries with experiments presents a 
stringent consistcincy test of the first order transition sce- 
nario. The first-order scenario was invoked by Ref. 43 
to explain the strongly enhanced longitudinal Coulomb 
drag for intermediate d/l, and it also has some support 
from exact-diagonalization study Note that we will 
only consider the case of negligible interlayer tunneling. 

The Clausius-Clapeyron relations are the results of 
matching the free energies of the two phases along the 
phase boundary. To be more specific, we denote the free 
energy density of the coherent and the incoherent phases 
to be Ec and Ei, and define 

f{5, Btou An, T) = Ec{5, Btot, An, T) - E^{6, Btot, An, T), 

where 5 = d/l, Btot is the total magnetic field coupled to 
electrons' physical spin, An = (m — n2)/2 is the density 
imbalance, T is the temperature. At any point along the 
phase boundary, we must have 



f{5c,Btot,An,T)^Q. 



(1) 



This equation can be viewed as defining the 5c at which 
the transition occurs. When one changes the total field 
by ABtot, the critical Sc{Btot, An, T) also changes by dSc 
when the filling factor is kept fixed at ftot = 1- Their 
relation is determined by 
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therefore the slope of the phase boundary is determined 
by the following ODE: 
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A crucial assumption of our work is that 
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where e^/(e/^) not only gives the correct units, but is 
the only energy scale that exists in this problem if we 
neglect the Landau Level mixing. 77 is a universal positive 
dimensionless constant. It is positive because / should 
be an increasing function of (5 = d/l, since the incoherent 
phase should be more and more energetically favorable 
with increasing d/l. In general, 77 could be a function of 
5 = d/l, i.e., rj{5) « //(Jq) + 0[(5 — 5q)/5q\, but since in 
experiments 5 docs not change much (ranging from 1.7 
to 2), {5 — 5q)/5o -C 1, we will assume r/ to be a constant 
for simplicity. 

Similar analysis also applies to finite temperature tran- 
sitions: 
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For density imbalance experiments, we will focus on the 
phase boundary near An = 0. First, note that by sym- 
metry 
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Thus, we need to expand / to second order in An: 
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and therefore 
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The above equations constitute the Clausius- 
Clapeyron relations for the bilayer quantum Hall 
systems. In the following sections, we will investigate 
whether the phase boundary shapes implied by Claiisius- 
Clapeyron relations are consistent with experiments, 
and whether a single universal parameter r] can explain 
all available experimental results. To obtain the detailed 
forms of free energy of both phases, we will primarily 
work with the pseiidospin ferromagnet description for 
the coherent quantum Hall phase and the Chern- Simons 
approach for the incoherent composite Fermi liquid 
phase. Spin transitions, finite temperature transitions, 
and density imbalance experiments are studied in Sec. 
II, III, and IV, respectively. Finally, we summarize and 
discuss our results in Sec. V. Some theoretical details 
are relegated to Appendices. 



II. SPIN TRANSITION EXPERIMENTS 

Ref. 26 and Ref. 30 have studied the effect of 
NMR/heat pulse and parallel magnetic field on the tran- 
sition tuned by d/l, respectively. In the experiment of 
Ref. 30, since the interlayer tunneling is negligible, the 
main effect of the parallel field is on the spins of electrons. 
Similarly, in the experiment of Ref. 26, NMR/heat pulse 
acts to depolarize the nuclei and therefore also changes 
the Zeeman field on the electrons through the hyperfine 
coupling. Thus, these two experiments can be analyzed 
in a similar fashion. Since we assume the coherent phase 
is spin polarized, the spin part of the coherent phase free 
energy is simply the Zeeman energy: 
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(9) 



where Nt is the total electron density of the two layers, 
B± is the perpendicular magnetic field, Btot is the total 
magnetic field coupled to electron spin, g = —0.44 is the 
g'-factor of the GaAs two dimensional electron gas, and 
fXB is the Bohr magneton. 

For the partially spin-polarized incoherent phase, the 
single layer free energy is 



E^ 
2 



1 

2X 



M2 - MBtot, 



(10) 



3 



where the magnetization 



M 



-j\9\l^B{n^ - nx) = \g\iJLB^n, (11) 



and X is the single layer spin susceptibility. The steady 
state is obtained by minimizing Ei with respect to M: 
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where the maximum magnetization Mmax and the field 
for full polarization Btot,p are given by 
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Plugging these forms of free energy into (3), we obtain 
an equation 
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Btot < Btot,p 
Btot > Btot,p 



(15) 



Note that the RHS also depends on Sc through B± which 
determines £. Eqn. (15) can be solved numerically to 
yield the 6c — Btot curve. For typical experimental pa- 
rameters, dSc/dBioi starts out to be positive when Btot 
is small, and continuously decreases to zero when 
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this is nothing but Eqn. (14) which determines the mag- 
netic field at which all composite fermions get polarized. 

It remains to determine the value of the composite 
fermion spin susceptibility x- This can be done if Btot 
and B± at which full polarization occurs are known, be- 
cause from Eqn. (14) or (16) we have 



X 



|g|MB-B-L,y 

4a 



tot,p<PO 



(17) 



where the subscript p denotes the point of full polariza- 
tion. In experimental and exact-diagonalization studies, 
one often parametrize x with the form of non-interacting 
Fermi gas with a "polarization mass" mj,^°'^^: 



X = 



(18) 



In the lowest-Landau-level approximation, ^ is the only 
relevant energy scale, and thus 



Pm„ ^ el 



(19) 
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FIG. 1: Total electron density deduced from the critical d/l = 
5c vs. the total magnetic field for the parallel magnetic field 
experiments. Open and solid circles arc experimental results 
(14) of Giudici et al}^° (c.f. FIG. 4a there). Solid line is our 
theoretical calculation with the fitting parameter rj = 0.8 • 
10~^. The boundary condition in our calculation is chosen as 
= 11 • 10i°cm-2 when B = lOT. 
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Therefore, presumably nip scales as 



nip = xnie 1 



(20) 



where nig is the vacuum electron mass, a; is a dimen- 
sionless number, B± is in units of Tesla. It is worth 
noting that unlike free electrons spin-susceptibility which 
is proportional to 1/nie, the susceptibility of composite 
fermions is proportional to nip and therefore to VB. The 
reason for this is that the Bohr magneton /i b depends on 
the bare mass of the electron, and therefore does not over- 
turn the proportionality to effective mass in the density 
of states factor of the susceptibility. 

For the parallel field experiment of Ref. 30, com- 
posite fermions get polarized at total density ntot = 
11- 10^°cm~^, tilting angle = 58°, which corresponds to 
Btot,p = 8.60T, Sx,p = 4.56T, x = 0.56 if we parametrize 
X in terms of the polarization mass mp. Then we solve 
the ODE (15) with the boundary condition at the high 
field endpoint {Btot = lOT, ntot = H • lO^^cm-^), and 
plot the ntot deduced from Sc vs. Btot in FIG .1. To time 
the result to resemble the experimental results in FIG. 
4a of Ref. 30, we get 



r/ = (0.8 ± 0.2) • 10" 



(21) 



where the error mainly comes from fitting errors, mean- 
ing a finite range of 77's make the 5c — Btot curve resemble 
the experimental result. 

For the NMR and heat pulse experiments of Ref. 26, 
the phase boundary before any perturbation is Sco ~ 
1.967, which correspond to B± = 3.26T. Ref. 26 
has estimated the effective nuclear magnetic field to be 
Bpf = — 0.17T, therefore the total effective magnetic 
field felt by electronic spin is Btot = B± + Bn- After 
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a heat pulse, nuclear spins arc depolarized, and is 
set to zero. Bfot is strengthened to _Bj_, and the phase 
boundary changes to 5c = 1.983. We can not determine 
the spin susceptibility or the polarization mass directly 
from experimental information, and therefore we use the 
numerical and experimental results from the literature 
rup = (0.7± 0.2)rne/BI with B in units of Tesla^°'52-55_ 
In this way, we obtain 



??« (1.3 ±0.4) • 10" 



(22) 



where the error mainly comes from uncertainty in the 
value of the polarization mass rUp. 

Note that our calculations in this section do not rely 
on the Chern-Simons description of composite fermions. 
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III. FINITE TEMPERATURE TRANSITION 
EXPERIMENTS 

Ref. 31 has studied the changes in critical Sc = d/l 
as a function of the temperature T. They found that 
the phase boundary moves to smaller d/l with higher T. 
When analyzing the temperature dependence of the tran- 
sition, one needs to include the entropy contributions to 
the free energy associated with various low energy excita- 
tions for both phases. In the interlayer-coherent quantum 
Hall phase, the only gapless excitation is the linearly dis- 
persing Goldstone mode, which corresponds to in-plane 
spin wave in the psciidospin language. Therefore, this 
mode dominates the temperature dependence of the free 
energy of the coherent phase (see 59). Denoting its ve- 
locity to be w, we have the free energy 



Ec{T) = ^rin(l - 6-'*"'=/'^) 
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We use the experimental result of Ref. 5 to estimate the 
value of V (which we assume to be a constant independent 
of S): 



1.4- lO^m- s"^ 



(25) 



For the incoherent phase, working in the Chern- 
Simons framework, we have contributions from compos- 
ite fermions as well as Chern-Simons gauge fields. The 
free energy is 



Ei = -TlnZ, 



(26) 



where the partition function Z contains both compos- 
ite fermion fields and Chern-Simons gauge fields of the 
two layers. Integrating out the composite fermions, we 
obtain^-^'^^ (see Appendix A for details) 



FIG. 2: The phase boundary d/l vs. the temperature T (in 
Kelvin) for the finite temperature experiments. Circles are ex- 
perimental results of Champag no et al.^^ (c.f. FIG. 2c there). 
Solid line is our theoretical calculation with the fitting param- 
eter ri = 0.7- 10^'*. The boundary condition in our calculation 
is chosen as 5c = 1.83 when T = 50mK. 



where Zq is the partition function for free fermions, and 

Z± = J Da±e--/''^^'*''^("±^^'"±/^\ (28) 

where a± are the in-phase and out-of-phase combinations 
of Chern-Simons gauge fields of the two layers, and the 
polarizations D^^ in the Coulomb gauge have the follow- 
ing form 
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where the index and 1 denote time and transverse com- 
ponent, respectively. 
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is linear combinations of intralayer and interlayer 
Coulomb interactions, F{q) is the finite thickness form 
factor^'"'^^, and IIqo and 11°^ are the fermion density and 
transverse current correlations functions, respectively: 
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m* is the activation mass of the composite fermions, and, 
as wc discuss below is different from the polarization mass 
rUp used in the previous section. Continuing the deriva- 
tion, 



Z = ZqZj^Z- 



(27) 
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(32) 
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where the free fermion part gives 
dEi 

^fermion d{TlnZo) _ 
dT dT 

and the gauge field parts give^^'^° 



(33) 
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A straightforward calculation following Rcf. 11 shows 
that in the zero-thickness approximation (form factor 
F{q) set to 1), 
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n is the single layer density of composite fermions, and 

kp = \/27rn. 

Finite thickness corrections to the form of Coulomb in- 
teraction is found to have negligible effect on the value 
of r], partly because it only affects the gauge field contri- 
bution which is itself dominated by the free composite- 
fermion-quasiparticle contribution for experimentally rel- 
evant temperatures and for the choice of m* discussed 
below. 

The value of the composite fermion mass m* is believed 
to be close to the value determined by the activation 
gaps of fractional quantum Hall phases away from v = 
-|^^2ii'i3,i5,56 Therefore, we use the experimental value of 
this activation mass determined from gap measurements 
in Refs. 61,62, which is 



0.2 ±0.02. 



(36) 



Note that in numerical calculations the activation mass is 
typically smaller than experimental value by about a fac- 
tor of 2ii'^-^'^^, but it is believed that the theoretical value 
should approach experimental value once finite thickness 
efifect, disorder and Landau level mixing are taken into 
account^"*"^^ . Therefore, we feel the use of experimental 
value stated above is more appropriate. Also note that 
the polarization mass nip we used in the previous section 
is diflFerent from the mass we use here. Conceptually, 
within the Landau Fermi liquid theory, the two masses 
are related by rup = m*/(l + i^Q ), Fq being the zeroth 
spin-asymmetric Landau parameter. 

Using this value of the mass along with the forms of free 
energy in Clausius-Clapeyron equation (5) , we get an 
ODE, which can be solved with the boundary condition 
that 5c = 1.83 when T 5GmK to yield the Sc - T 
curve plotted in FIG. 2. To make this curve resemble the 
experimental result of Ref. 31, we have set 



77= (0.7 ±0.2) • 10" 



(37) 



where the error mainly comes from the uncertainty in 
the value of the activation mass m* and also the fitting 
error, meaning a finite range of 77 's make the 5c — T curve 
resemble the experimental result. 

In the above calculation, we assumed that the com- 
posite Fermi liquid is spin-unpolarized, and one might 
wonder how partial spin-polarization would afi'ect the re- 
sult. Because the free fermion contribution dominates 
dEi /dT and it is proportional to the density of states of 
composite fermions, our results would stay the same for 
partially-polarized composite Fermi liquid. 



IV. DENSITY IMBALANCE EXPERIMENTS 

Refs. 24,32 have studied the dependence of the criti- 
cal 5c = d/l on the density imbalance between the layers. 
They observed that at small imbalance, the phase bound- 
ary has a quadratic dependence on the density imbalance, 
and the coherent quantum Hall phase survives at higher 
d/l with larger imbalance. 

Denoting the density of the two layers ni,2, a density 
imbalance between the layers. 



An = 



ni - n2 



(38) 



costs an energy which includes a dominating geometrical 
capacitance term and quantum mechanical corrections. 

This is true for both phases. For the coherent phase, we 
follow Ref. 2 to obtain the free energy density to be 
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where l-ne^d/ e is the geometrical capacitance term, while 
/3m, E is the exchange contribution which tends to off- 
set the geometrical capacitance term. Here, F^(g') = 
V{q) ~ U{q), V{q) = ^F{q) is the intralayer Coulomb 

interaction, F{q) is the finite thickness form factor^'''^^, 
U{q) = V{q)eri'^ is the interlayer Coulomb interaction, 
and h{q) = —2TTpcxp{—q^P/2) is the pair distribution 
function of the Halperin (1,1,1) wavefunction. 

The free energy density of the incoherent phase is (see 
Appendix B for details) 
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where 



^=1U I™ lim(Pl,9,a;Pl,-g-,-a;), 
pA g-s-O aj->0 

k' = ^ lini lim (pi,,-,„p2 -g-a;), 



(41) 



where j3 is the inverse of temperature, A is the area of 
the sample, pi^2 are the composite fermion density of each 
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layer. Treating the Coulomb interaction within RPA, we 
obtain (see Appendix B for details) 
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where n is the w — 0,(7 — >■ limit of the 1-particle- 
irreducible density response function, namely compress- 
ibility, of a single-layer composite Fermi liquid. Plugging 
(42) into (40), one obtains the energy cost of uniform 
density imbalance in the incoherent phase: 



1 2TTe^d 
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From the Clausius-Clapcyron equation (8), the geomet- 
rical capacitance term of the two phases cancels out, and 
we have 



(44) 



Since k is the single layer compressibility, it is con- 
nected to the ground state energy per area of the com- 
posite Fermi liquid Ecs via 



-1 _ d^Eps 



(45) 



Note that our definition of the compressibility is slightly 
different from some literature where = 'n? ^g^f ^ are 
used instead. 

Alternatively, treating the Chern-Simons interaction 
within RPA (see Appendix B for details), we obtain 



where 



Kq = 



(46) 



(47) 



is the compressibility without the Chcrn-Simons inter- 
action, Fq is the zeroth Landau parameter in the spin- 
symmetric channel, and 



Xd 



127rm* 

is the Landau diamagnetic susceptibility. Therefore 



(48) 
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Clearly, we can identify the four terms as free fermion 
contribution, exchange/correlation effect, Landau dia- 
magnetism for Chern-Simons flux^^, and geometric ca- 
pacitance term, respectively. 

Although the Chern-Simons expression of k Eqn. (49) 
offers valuable physical insight into its structure, the pre- 
cise value of the parameters m^, Xd, and especially Fq 
are not very well understood. The best way to estimate k 



is to use its connection with ground state energy density 
Eqs of composite Fermi liquid (45). In the zero-thickness 
approximation. Park et al.^^ have estimated the vahic of 
Egs for spin unpolarized composite Fermi liquid to be 



thus 



Egs = -0.4695— n, 
el 



-0.4695 • 37r—Z, 
e 



(50) 



(51) 



where n is the single layer density of composite fermions, 

and / is the magnetic length. 

Using this value of k^^ and the zero-thickness form of 
Coulomb interaction to calculate the coherent phase ex- 
change term l3m,E (because the numerical result for Egs 
of the incoherent phase quoted above from Ref. 68 was 
also done with zero thickness), and extracting the curva- 
ture ^^^°2) from experiments, we readily obtain the value 
of T]. This result does not depend on the Chern-Simons 
description of composite fermions. We have plotted in 
FIG. 3 the values of rj extracted from density imbalance 
experiments as well as those determined from spin transi- 
tion and finite temperature transition experiments. The 
error bars for the density imbalance experiments mainly 
come from fitting errors. 

Note that the main effect of the finite thickness cor- 
rection to the form of Coulomb interaction is to reduce 
the exchange terms of both phases. Since the value of 
T] is related to the difference between the exchange term 
of the two phases, we do not expect the result of rj to 
sensitively depend on this effect. Nevertheless, we can 
include it in the Chern-Simons treatment of k. We use 
the activation mass m* = 0.2me\/B± estimated in Sec. 
Ill as the value of m*, set Xd = —l/{12^^m.^,), and use the 
Hubbard approximation to estimate Fq . In the Hubbard 
approximation, the exchange effect is taken into account 
by introducing a many-body local field factor G{q) = 

«/(2y?TfcF). and FS = -^lim,^oV{q)G{q). Thus, 
we obtain from Eqn. (46) 
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Using this value of and the finite-thickness form 

of Coulomb interaction to calculate the coherent phase 
exchange term (3m, e, we have calculated the values of rj 
from density imbalance experiments which turned out to 
be extremely close to the results obtained earlier in FIG. 
3. 

Comments about the value of the compressibility in 
the composite Fermi liquid phase are in order. First, In 
Ref. 69, the compressibility of a single layer 2DEG at 
zero field was studied in detail, and it was found that 
aside from the well-known density-of-states contribution 
and exchange contribution to the compressibility, there 
is a third contribution coming from the so-called Hartree 
band-bending effect due to the influence of the finite 
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FIG. 3: Summary of the value of ij extracted from various ex- 
periments. Experiment 1: parallel field experiment of Ref. 30. 

Experiment 2: NMR/hcat pulse experiment of Ref. 26. Ex- 
periment 3: finite temperature transition experiment of Ref. 

31. Experiment 4 to 7: density imbalance experiments of Ref. 

32, with T = 55mK, 85mK, 125mK, 200mK. Experiment 8 
and 9: density imbalance experiments of Ref. 24 with phase 
boundary determined by Hall drag and tunneling. To obtain 
this result we used the rmmcrical result of Ref. 68 for unpolar- 
izcd composite Fermi liquid ground state energy to estimate 

. The horizontal line is the average value of rj weighted 
by inverse of error square, which is ~ (1 ± 0.1) x 10~^. 



quantum well width on the out-of-plane direction of elec- 
tron wavcfunction. For the bilaycr system studied here, 
we expect a similar effect on the composite Fermi liquid 
compressibility in the incoherent phase and on Pm,E 
for the coherent quantum Hall phase as well. A quanti- 
tative analysis of this effect and its impact on the density 
imbalance experiments is beyond the scope of this paper, 
and we simply note that the Hartree band-bending ef- 
fect is essentially a single-particle effect^'', and therefore 
it will contribute equally to and l3m,E- To obtain the 
value of r] from Eqn. (44), we only need the difference 
between and /3m, b, and therefore we do not expect 
the Hartree band-bending effect to modify our results. 
Second, quenched disorder acts to broaden the Landau 
levels and therefore adds a positive contribution to the 
compressibility. This could account for the close-to-zero 
compressibility measured by Ref. 69. Again, this effect 
is likely to be similar for both phases, and we do not 
expect disorder to affect the difference between and 
/3m, E appreciably. Nevertheless, disorder is important in 
smearing the first order transition into a continuous one 
(see discussion in Sec. V). 

We assumed that the composite Fermi liquid is unpo- 
larized above, but again we do not expect partial polar- 
ization to affect our results strongly. For (51), Park et 
al.^^ also reported the ground state energy for polarized 
composite Fermi liquid to be very close to the unpolar- 



ized one quoted above: 

g2 

Epolarized = -0.4656— n, (53) 

and therefore our results would also stay very close. 
In the Chern-Simons treatment (49) and(52), since the 
Chern-Simons fields couple to both spins and the den- 
sity and current response function stays the same for 
partially-polarized and unpolarized composite Fermi liq- 
uids, our calculation also remains valid (see Appendix 
B). 



V. SUMMARY AND DISCUSSION 

To summarize, we derived the Clausius-Clapeyron re- 
lations [Eqn. (3, 5, 8)] for the phase transition tuned by 
d/l in bilayer vtot = 1 quantum Hall system, assuming 
that it is a first-order transition between spin-polarized 
coherent quantum Hall state and spin partially-polarized 
composite-fermion Fermi liquid state. In Sec. II, we stud- 
ied the changes of phase boundary {d/l)c when the mag- 
netic field coupled to spin is changed by either NMR/heat 
pulse or parallel magnetic field. The phase boundary as a 
function of temperature was studied in Sec. HI. The tem- 
perature dependence of free energy in the coherent quan- 
tum Hall phase is dominated by the linearly-dispersing 
Goldstone mode, while the incoherent composite Fermi 
liquid phase has contributions from both fermions and 
gauge fields. In Sec. IV, we investigated the changes of 
phase boundary when there is density imbalance between 
the two layers. We use the result of Ref. 2 for the free 
energy cost of density imbalance in the coherent quan- 
tum Hall phase. The free energy for the incoherent phase 
is shown to be connected to the compressibility of single 
layer composite Fermi liquid. 

Our main goal was to check the consistency of the 
Clausius-Clapeyron relation with the observed transition. 
Each experiment which observes the change in {d/£)c due 
to changing another parameter in the system indicates a 
value for r], as defined in Eq. (4); all values should agree. 

In FIG. 3, we have plotted the values of t; determined 
from spin transition, finite temperature transition, and 
density imbalance transition experiments. The horizon- 
tal line is the average value of 77 weighted by inverse 
of error square, i.e., the maximum likelihood estima- 
tor of T]. One can see that, indeed, all nine values of 
Tj extracted from various experiments roughly lie in the 
range 1 ~ 2 x 10~^, and the weighted average value of 
?7 = (1 ± 0.1) • 10~^ is roughly within all the error bars. 
Our analysis, therefore, confirms the consistency for the 
scenario of a direct first-order phase transition between 
coherent quantum-Hall phase and incoherent composite 
Fermi-liquid phase. Furthermore, the analysis provides a 
unified framework within which we can understand the 
observed phase boundaries for several distinct experi- 
ments. 
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FIG. 4: (Color online.) Standard deviation of rj among 
various experiments divided by their average value (which 
measures the goodness of agreement between rjs) within the 
Chern-Simons framework as a function of the composite 
fermion mass and the exchange contribution to [see Eq. 
(52)]. Horizontal axis: composite fermion mass in units of 
me\/B±, me being the vacuum electron mass, B± is in units 
of Tesla. Vertical axis: exchange contribution to which 
is FQ-K/rrit (in units of e^l/e, I being the magnetic length). 
Grey color denotes the region where at least one of the ry's be- 
comes negative, thus unphysical. The horizontal line denotes 
the Hubbard approximation to the exchange effect (— \/27r). 
The vertical line denotes the experimental value of the activa- 
tion mass m, ~ Q.2m^^ B± , which is the value of composite 
fermion mass we used in calculations for FIG. 2 and FIG. 3. 



In Sec. IV, we also worked in the Chern-Simons de- 
scription of composite fermions [i.e. Eq. (52)] in addi- 
tion to our treatment [i.e. Eq. (51)] using the numerical 
results of Ref. 68, and we obtained extremely similar 
results. Stepping back a little from that analysis with 
the Chern-Simons treatment, one can pretend ignorance 
of any knowledge of the parameters including the effec- 
tive mass m* and the exchange contribution to k^^, and 
ask what values of them would give good agreement be- 
tween the values of rj extracted from experiments. We 
have plotted the standard deviation of rj extracted from 
various experiments divided by the their average value in 
FIG. 4 as a function of the composite fermion mass (in 
units of me\/B±) and the exchange contribution to k^^, 
which is i^Qvr/m, (in units of e^l/e). The finite thick- 
ness form of the Coulomb interaction is used in calculat- 
ing the coherent-phase exchange term when producing 
this plot. Grey color denotes the region where at least 
one of the ?7's becomes negative, thus unphysical, while 
dark blue denotes parameter regimes which give rise to 
good agreement among yy's extracted from different ex- 
periments. The horizontal line denotes the Hubbard ap- 
proximation to the exchange effect (\/27r), while he verti- 



cal line denotes the experimental value of the activation 
mass TO* 0.2TOe-\/i3j_. 

We have not explicitly discuss the role of disor- 
der, which is always present in the samples. Disor- 
der will bring spatial fluctuations into some variables in 
the Clausius-Clapeyron equations we have derived, and 
therefore smear the flrst order transition into a continu- 
ous one, as observed in experiments. Roughly speaking, 
the analysis we have performed in this work applies to 
the spatially averaged quantities. For example, with dis- 
order, the RHS of the Clausius-Clapeyron equation for 
spin transitions (15) will acquire spatial dependence most 
likely through a spatially fluctuating spin susceptibility 
X- 



Thus, one can take the spatial average of both sides and 
study how the averaged critical 5c changes with Btot, as 
we did in this work. Furthermore, one can also take the 
standard deviation of both sides of (54), and conclude 
that the width of the phase transition, which is the stan- 
dard deviation of 5c, grows with Btot assuming the stan- 
dard deviation of x(a;) does not change appreciably with 
Btot- One can also study the finite temperature tran- 
sition in a similar way. Because there the free fermion 
term (33) dominates, one can conclude that the transi- 
tion becomes wider at higher temperature, if one assumes 
the composite fermion mass to* has some temperature- 
independent spatial variation. This is in accord with the 
experimental observation of Ref. 31. 

A major question which is not directly addressed in 
our analysis is the possibility of a continuous quantum 
crossover between the coherent and incoherent phases 
(see, e.g., Refs. 47,48). If indeed no real thermodynamic 
singularity exists even in the clean case, then there is 
no reason for the Clausius-Clapeyron relations to hold as 
well as we find they do. Nonetheless, there is also no con- 
tradiction in them holding where no first-order transition 
exists. In this case, however, we can draw the conclusion 
that the crossover region between the two phases must 
be very narrow, such that it approximates a smeared 
thermodynamic singularity (just as disorder would widen 
a thermodynamic singularity) and therefore follows the 
Clausius-Clapeyron relations we presented here for the 
unmixed phases. In other words, a good agreement with 
the relations indicates that already at regions in param- 
eter space close to the transition, the thermodynamic 
functions of the pure coherent and pure incoherent phases 
apply, and they indicate a smeared phase transition line. 

Additional outstanding questions which we did not ad- 
dress, but are noteworthy are as follows. First, a ther- 
modynamic phase transition between the coherent and 
incoherent phases does not have to be first order at high 
Zeeman fields when both phases are spin-polarized; a 
second-order phase transition is not ruled out a priori. 
Future experiments should clarify this issue (see the re- 
cent experiments of Refs. 34,35). In addition, for the 
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density imbalance transitions, we have mainly focused 
on the regime of small imbalance, while the experiments 
of Ref. 32 have studied the case of large imbalance, e.g., 
Au = v\ — V2 < 0.4. The interlayer incoherent phase 
in that regime could be two decoupled single-layer frac- 
tional quantum Hall phase. It would be very interesting 
to see if a similar Clausius-Clapeyron equation can de- 
scribe the phase transition in that case. Finally, although 
our assumption (4) is very natural on qualitative ground, 
a microscopic derivation of this quantity would be very 
useful. 
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Appendix A: DERIVATION OF THE 
TEMPERATURE DEPENDENCE OF THE 
COMPOSITE FERMI LIQUID FREE ENERGY 

Within the Ghern- Simons description of the 
composite- fermion Fermi liquid at v = 1/2, we 
have the following partition function of the system: 

Z = JvaiVa2ViPi„Vi;2ae-^'^^'^'^'^, (Al) 

where 



^=l,2 



+ Ti—i'L - a„) Ipr, 



2m 

+^ j d2a;VLWV'n.(a;)^(a;-a;')VL'(a=')V'na'(^')} 

+ j d^x''4)\„{x)'ll}ia{x)U{x - x')-li)\^,{x')-li)2a'{x'), 

(A2) 

where tpna is the composite fermion fields in the n'th 
layer with spin a, V and U are the intralayer and inter- 
layer Coulomb interaction, respectively. Here, are 
the fluctuations of the Ghern-Simons gauge fields in the 
n'th layer from its saddle point value which cancels the 
external magnetic field exactly, and fjb = 0,1,2 are the 
time and two spatial coordinates, respectively. Integrat- 
ing out o„,0i one obtains the expected constraints 



Following Ref. 11, we make use of this constraint and 
replace ipj^^ipna in Coulomb interaction terms by V x 
a„/(47r). Next, we define 



^ 2 ' 

and reorganize C as 
^ = + ^cs, 

r /t (a .a+o-|-a_o\ , 



(A4) 



1<T 



^ HV-p. + .--)/2)% ,^ 

Zm 

,t (-iW - (d+-a-)/2f , 

Ccs = --^a+.e^'^^d^a+x - -^-a.^e^'^^d^a-x 
loTT Idtt 

+ j d'^'[V X a+{x)]V+{x - x')[V x a+{x')] 

^ ^ j d^x'[V X a-{x)]V-{x - a;')[V x a-{x% 

(A5) 



+ 



2 (47r)2 

Denoting the free fermion partition function to be 

Zo = J I'ViaI'V2a exp J dTd'^xCf{a± = 0)^ , 

(A6) 

and following standard methods^'^ to integrate out com- 
posite fermion fields ipna, we obtain 

Z = ZoZ+Z-, (A7) 

where Zq is the partition function for free fermions, and 

Z+ = /r>a+e-i''*^'^'^("±^^'"±/2), (A8) 



In Coulomb gauge, one can treat the polarizations 
as 2 X 2 matrices, with index and 1 to be the time and 
transverse component, respectively. Thus, D^^ take the 
following form: 



D 



-1 ifK 



00 47r 

± -2l ^ no 



(A9) 



V X a„ = Anipl^'tpna. 



(A3) 



where Hqq and H^^ are the density and transverse cur- 
rent correlation functions of free fermions resulted from 
integrating out composite fermion fields. Thus, the free 
energy is given by 

Ei = -TlnZ = -TlnZo-TlnZ+-TlnZ_, (AlO) 

and the rest of the steps are given in Section HI. 
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Appendix B: DERIVATION OF THE FREE 
ENERGY FOR DENSITY IMBALANCE IN 
COMPOSITE FERMI LIQUID PHASE 

Starting from action (Al) or any other action for com- 
posite fermions, we integrate out all fluctuating flelds and 
obtain 



Z = exp 



ui9l,q,u'Pl,-q,-w 



(Bl) 



where 

= -^{Pl,q,U)Pl,-q,-U}) = -p^{P2,q,uP2,-q,-u) , 
K'q,ui = ■^{Pl,q,U}P2,-q-u}), 

(B2) 

Pj is the composite fermion density of the j'th layer, ^ 
is the inverse of the temperature, A is the area of the 
sample, and <pj,q,ui is the Fourier-transformed potential 
in the j'th layer. For a constant potential (f)j {j = 1,2), 
we have 



and the grand potential O is 

A 



(B3) 



n = -T In Z = - ^ + iCcPl + 2k'(j)^(j>2) , (B4) 



where 



K = lim lim K^,,,, K' = lim lim K';^, ,. 



(B5) 



The density in each layer is 

1 dfl 
Ad4>2 

Finally, the free energy is obtained via a Legendre trans- 
formation 



(B6) 



F = 0. + (pmiA + <p2n2A 

A f {m - 722)^ (ni + 712) 



+ 



4:\K-K' K + K 



(B7) 



Within the RPA treatment of the Coulomb interaction, 
the full density response function K is related to its one- 
particle-irreducible (IPI) counterpart 11 (which neglects 
the long range Coulomb interaction) by 



K-'^ = U-^ + V, 



(B8) 



where K, 11, and V are 2x2 matrices in the layer-index 
space: 



V 



V U 
U V 



Hoo 
Hoo 



(B9) 



Here, V and U are intralayer and interlaycr Coulomb 
interaction potential, respectively, and lloo in the static 
uniform limit gives the single layer compressibility k: 



K = lim lim IIoo. 



(BIO) 



Solving (B8), we have 



22 



noo(i + nooi^) 



K 



12 



(i + nooV^)2-n2o[/2- 

Given the form of Coulomb interactions 

27re2 



-F{q\ U{q) = V{q)e-'^\ (Bll) 



V{q) 



and the fact that the finite thickness form factor F{q) — >■ 
1 as g — >■ 0, in the limit w — )■ and g — >■ 0, the denomi- 
nators of Kii, K22, and K12 become 



{1 + IlooVf-iIloofU^ 

Aire'^K ( 27re^K(i^ 
I 1 -I I , as w 



Therefore in this limit 

K = lim lim Kn = 

k' = lim Um K12 = 



0,q->-0. 



(B12) 



2 (1 -h 2Tre^Kd/e) 

K 



(B13) 



q^o " 2 (1 -I- ^ne^Kd/e) ' 
and the imbalance part of the free energy density is 



Ei = lim lim — 



g-^O w^o K — K' 

1 27re2d\ , 2 

- + An^, 

K el 



(B14) 



as shown in Section IV. This result does not depend 

on the Chcrn-Simons description of composite fermions. 
Note also that the total compressibility K + K' vanishes 
linearly mqasq^Q due to the long-range nature of the 
Coulomb interaction, similar to the single layer case as 
analyzed by Halperin et al}^ . 

To calculate the single layer compressibility k within 
the Chern-Simons framework, we have the following RPA 
equation: 



(n) 



(n") 



(B15) 



where C is the propagator of the Cherns-Simons field, 
and n*^ is the correlation functions without the Chern- 
Simons interaction. Wc work in the Coulomb gauge and 
treat 11, 11°, and C as 2 x 2 matrices in the space of 
density and transverse current. In the static and long 
wavelength limit, we have 



no 



KO 

Xdq^ 



C= — 

q 



i 
-i 



(B16) 
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where kq = TO*/[7r7i^(l + Fq)] is the density response 
function neglecting Chern-Simons interaction, and Xd is 
the Landau diamagnetic susceptibiUty. Hence, 

= Ko ' - IfiTT^Xd, (B17) 

as shown in Section. IV. Note that these results are the 
same for unpolarized and partially-polarized composite 



Fermi liquids, because (B15) is valid in any case since 
Chern-Simons fields couple to both spins, and the value of 
kq and Xd in (B16) stays the same for partially-polarized 
composite Fermi liquid. The value of Fq in the Hub- 
bard approximation treatment is also roughly the same 
for partially-polarized and unpolarized composite Fermi 
liquids. 
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